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Abstract. Simple spherical, non-rotating stellar models are inadequate when describing real 
stars in the limit of very fast rotation: Both the observable spectrum and the geometrical 
shape of the star deviate strongly from simple models. We attempt to approach the problem 
of modeling geometrically distorted, rapidly rotating stars from a new angle: By constructing 
distorted geometrical models and integrating standard stellar models with varying temperature, 
gravity, and abundances, over the entire surface, we attempt a semi-empirical approach to 
modeling. Here we present our methodology, and present simple examples of applications. 



1. Introduction 

For many years, stellar atmospheric models and spectrum synthesis has been synonymous with 
the use of standard ID atmospheric models. Lately, 3D models have added another dimension 
of realism to this picture. However, all models in use today are local models, treating only a 
small part of the stellar surface. Calculation of synthetic spectra then relies on assuming that 
the stellar atmosphere looks the same everywhere. Real stellar atmospheres are of course not 
the same everywhere. In this contribution we wish to highlight an aspect of stellar spectral 
synthesis that is mostly overlooked, namely the deviations from spherical symmetry. 

The problems associated with stars that are not spherically symmetric were first addressed 
by von Zeipel [1; 2], who showed that stars that were flattened due to fast rotation would 
have a latitude dependent brightness distribution. However, the observed spectrum of a star 
will be influenced by many different additional effects from differential rotation, turbulence, 
temperature and gravity variations — all latitude dependent. Not to mention effects like spots, 
plagues, etc. that are both latitude and longitude dependent. 

2. Methodology 

Ideally, one should be able to take an observed spectrum and perform an inverse analysis, that 
is, from the spectrum infer the exact distributions of temparature, gravity, abundances, etc., as 
well as the shape of the star. However, such capabilities are far into the future since we still 
don't understand fully all the effects going into the structure of a star and the formation of its 
spectrum. 



Figure 1. Example of a tiled spherical surface with 22 latitudinal bands (N = 11) before the 
flattening coordinate transformation is applied. Note that the shape of the triangular tiles varies 
slightly across the surface. 

With this work, we are attempting a different approach to learn more about these effects, 
that determine the structure of geometrically distorted stars and the possibilities of observing 
those effects. Our approach is: 

• We assume a shape of the star - any shape. 

• We assume (any) distributions of temperature, gravity, abudances, etc. across the surface. 

• We assume any rotation law on the form v{9). 

• At any point on the stellar surface we then compute a resulting intensity spectrum, in the 
direction of the observer. For this, any kind of local stellar model atmosphere can be used, 
whether ATLAS, MARCS, or state-of-the-art 3D models. 

• We then sum up all the contributions from the visible part of the star into an "observed" 
synthetic spectrum. 

This is an exploratory approach, aimed at providing more insight into the individual effects 
involed in shaping the spectra of real stars, rather that an attempt at producing realistic models 
of real stars, although that is of course the long term aim. 

We want to avoid dependence on analytical descriptions of the parameters, which means that 
we must be able to freely assign any value of any parameter to any point on the surface. For this 
reason we chose to create a tiled surface, i.e., we cover the entire surface with plane segments to 
which we can assign any combination of parameters we wish. In principle the only limit to the 
fine graininess of the surface is the speed of the computer hardware. Since a sphere cannot be 
tiled with equal-area regular polygons, we chose instead to employ triangular tiling where the 
individual tiles have slightly different surface areas. The tiles are arranged in latitudinal bands, 
where N is the number of bands per hemisphere. See Fig. 1 for an example. 
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Figure 2. Flow chart describing the process. The blue boxes represent the software modules 
described in this paper. The box SupeRotate is where the actual spectral synthesis is done. The 
database of synthetic intensity spectra can be constructed using any current model atmosphere 
grid. 

The full procedure of what we have named "SupeRotate" is shown schematically in Fig. 2. 
The details will be treated in a later paper. 

3. Results 

As a zero-order test of the code we first simulated a non-flattened star with no variation of 
parameters across the surface and compared this to a standard synthetic spectrum derived from 
the same set of models. We tested this for models with varying rotational velocities, assuming 
rigid rotation. The results were practically identical for N > 9. 

As more interesting tests we have simulated two cases which resemble well known stars, 
namely Altair and FK Com. However, these models are intended to test the model sanity rather 
than produce actual physical results on the two stars. 

3.1. Models of "Altair" 

Altair is an extensively studied flattened, fast rotating early type star. Its oblateness or 
flattening of r po \/r eq = 0.80 was first demonstrated by van Belle et al. [3] with subsequent 
studies confirming these findings [4; 5]. For this study, for the purpose of illustration, we are 
modeling Altair with an equatorial rotational velocity of 80 kms -1 as opposed to the true 
velocity of about 270 kms -1 . Thus, we refer to our models as "Altair" rather than Altair. 
We have calculated synthetic spectra for four different inclinations (see Fig. 3) with 
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Figure 3. Close up of a section of the resulting synthetic spectrum of "Altair" with T e ff and 
logg varying from pole to equator, and at four different values of the inclination; i = 0° (black), 
i = 30° (blue), i = 60° (magenta), i = 90° (red). 



Table 1. Parameters used in the simulations of "Altair". The model is then used to calculate 
synthetic spectra for four values of the inclination; i = 0°, 30°, 60°, 90°. 



Parameter 


Range (pole — equator) 


^eff 


7300 K — 6900 K 


log# 


4.1 — 3.9 



atmospheric parameters varying linearly in latitude from pole to equator, 
used are summarized in Table 1. 

Further tests on "Altair" is planned and will be reported elsewhere. 
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3.2. Models of "FK Com" 

FK Com is, with its usinz = 160 kms" 1 , the fastest rotating and most active of the FK Comae 
stars (see [6] and references therein). These are a small group of magnetically very active 
late-type giants with rotation periods of only a few days. The photometric and spectroscopic 
characteristics of FK Comae stars are very similar to those of the very active RS CVn stars, but 
the RS CVn stars are close binary systems while FK Comae stars are single, possibly the result 
of a merging of the components of a RS CVn system. 



Table 2. Parameters used in the simulations of "FK Com" as well as the average parameters 
from [7]. The spectra are all calculated for the known inclination i = 60°. 



Parameter Range (pole — equator) Standard model 



T eff 5500 K — 4600 K 

logg 3.8 — 2.8 
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Figure 4. Observed spectrum of FK Com (black) along with our synthetic spectrum from 
SupeRotate (orange) and from standard synthesis (blue). 



The parameters of FK Com is well known and lies within the model grid we have prepared 
for our testing phase. Though it is believed that temperature, gravity and indeed abundances 
are varying widely over the surface, we have employed parameters varying smoothly with 9 only. 
Thus, our model is "FK Com" rather than FK Com. Table 2 lists the parameters used, as well 
as the standard average stellar model atmosphere consensus for FK Com. In our calculation 
we used the known vsini and inclination i = 60° as well as an assumed flattening of 0.80. The 
resemblance with the observed spectrum is good as is evident from Fig. 4. 

In Fig. 5 we show the results of the flattened model campared with the standard model for 
a small spectral range. It is clear that the differences are quite small and possibly too small to 
be observed, in particular when one considers that the surface of FK Com is a mix of different 
temperatures, abundance patches, active regions, etc. However, what this code can provide is 
firm ideas about the relative significance of all these effects and what changes might be induced 
by variations of any of these parameters. 
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Figure 5. Synthetic spectrum from SupeRotate (orange) and from standard synthesis (blue) 
in a 70 A interval. 

4. Conclusion 

In this contribution we have explored a new modeling approach which is aimed primarily at 
gaining insight into the underlying individual effects that are shaping the spectra of stars with 
non-spherical and non-homogenious surfaces. 

Our models are currently at the "toy model" stage, suited for an exploratory approach. Our 
models are easy to use, flexible, and highly adaptable. Following further testing, the software 
will be made freely available to the community. Meanwhile, we welcome requests and suggestions 
for suitable projects which we will undertake as extended testing. 
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